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Abstract 

The production of heavy sterile neutrinos from 7r~, K~ decay at rest yields charged leptons with 
negative helicity {positive for tt + ,K + ). We obtain the branching ratio for this process and argue 
that a Stern-Gerlach filter with a magnetic field gradient leads to spatially separated domains of 
both helicity components with abundances determined by the branching ratio. Complemented 
with a search of the monochromatic peak, this setup can yield both the mass and mixing angles for 
sterile neutrinos with masses in the range 3 MeV < m s < 414 MeV in next generation high intensity 
experiments. We also study oscillations of light Dirac and Majorana sterile neutrinos with m s ~ eV 
produced in meson decays including decoherence aspects arising from lifetime effects of the decaying 
mesons and the stopping distance of the charged lepton in short baseline experiments. We obtain 
the transition probability from production to detection via charged current interactions including 
these decoherence effects for 3 + 1 and 3 + 2 scenarios, also studying |AL| = 2 transitions from 
Q_[ v •<-)■ v oscillations for Majorana neutrinos and the impact of these effects on the determination 

of CP-violating amplitudes. We argue that decoherence effects are important in current short 



baseline accelerator experiments, leading to an underestimate of masses, mixing and CP-violating 
angles. At MiniBooNE/SciBooNE we estimate that these effects lead to an ~ 15% underestimate 
for sterile neutrino masses m s > 3eV. We argue that reactor and current short baseline accelerator 
experiments are fundamentally different and suggest that in future high intensity experiments with 
neutrinos produced from n, K decay at rest, stopping the charged leptons on distances much smaller 
than the decay length of the parent meson suppresses considerably these decoherence effects. 
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I. INTRODUCTION 



Neutrino masses, mixing and oscillations are the clearest evidence yet of physics beyond 
the standard model [2-0] . They provide an explanation of the solar neutrino problem p-0] 
and have important phenomenological jH 0, B, i-El, astrophysical @, [H, Q and cosmo- 



logical |15| consequences. A remarkable series of experiments have confirmed mixing and 
oscillations among three "active" neutrinos with 5m 2 = 10~ 4 — 10 _3 eV 2 for atmospheric and 
solar oscillations respectively flfj] . 

However, several experimental hints have been accumulating that cannot be interpreted 
within the "standard paradigm" of mixing and oscillations among three "active" neutrinos 
with 5m 2 ~ 10~ 4 — 10 -3 . Early results from the LSND experiment [13] have recently been 
confirmed by MiniBooNE running in antineutrino mode [18] both suggesting the possibil- 



ity of new "sterile" 
collaboration! H 



neutrinos with 5m 2 ~ eV 



on the combined — > v e and 
neutrino oscillations with 0.01 < Am 2 < 1.0 eV 2 



The latest report from the MiniBooNE 
— > V e appearance data is consistent with 
This is consistent with the evidence from 
LSND antineutrino oscillations [13], which bolsters the case for the existence of sterile neu- 
trinos; however, combined MiniBooNE/SciBooNE analysis (20| of the disappearance data 
are consistent with no short baseline disappearance of V^. Recently, a re-examination of the 
antineutrino ffux(2l| in anticipation of the Double Chooz reactor experiment resulted in a 
small increase in the flux of about 3.5% for reactor experiments leading to a larger deficit 
of 5.7% suggesting a reactor anomaly^^. If this deficit is the result of neutrino mixing and 
oscillation with baselines L < 10 — 100 m, it requires the existence of at least one sterile neu- 
trino with 5m 2 > 1.5 eV 2 and mixing amplitude sin 2 (2#) ~ 0.115|22]. Taken together these 
results may be explained by models that incorporate one or more sterile neutrinos that mix 



with the active ones[23H30j including perhaps non-standard inter actions [3 11]; although, there 



is some tension in the sterile neutrino interpretation of short-baseline anomalies 32]. A com- 



prehensive review of short baseline oscillation experiments summarizes their interpretation 
in terms of one or more generations of sterile neutrinos 33 . 



Recently it has been pointed out that the presence of sterile neutrinos may induce a 
modification of the recently measured angle #13 34|, |35| . 



Hints for the existence of sterile neutrinos also emerge from cosmology. The analysis 
of the cosmic microwave background anisotropies by WMAP|36] suggests that the effective 
number of neutrino species is N e ff = 4.34 ± 0.86 and XX m ^) < 0.58 6^, bolstering the 
case for sterile neutrino (s) with m < eV. Complementary cosmological data suggests that 
N e ff > 3 at the 95% confidence level [3 71]; although, accommodating an eV sterile neutrino 
requires a reassessment of other cosmological parameters [38]. For recent reviews on "light" 
sterile neutrinos see ref. [39]. Furthermore, sterile neutrinos with masses in the ~ keV range 
may also be suitable warm dark matter candidates Mil. While this mass range is much 
larger than the hints emerging from accelerator and reactor experiments, it is intriguingly 
noteworthy that the ratio of this "heavier" species to the "lighter" ones (m ~ 1 eV) is of the 
same order as that in the mass hierarchy of charged leptons. 

When taken together, these emerging hints motivate several experimental proposals to 
search for sterile neutrinos (see the reviews in ref. [39]). Various experimental searches have 
been proposed, such as Higgs decay and matter interactions of relic sterile neutrinos [461]. the 
end point of /3-decay in 187 Re with a value of Q = 2.5 keV(47l. |48| (although the statistics 
will be hindered by the long lifetime of the source ~ 4.3 x 10 10 years), and electron capture 
decays of 163 Ho -> 163 Dy{49] with a Q-value ~ 2.2 keV — 2.8 keV. More recently, the focus 
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has turned on the possible new facilities at the "intensity frontier" such as project X at 
Fermilab[50], alternative high intensity sources |51| and recent proposals to study sterile- 



active oscillations with pion and kaon decay at rest (DAR) [52|, |53j or muons from a storage 



ring[54[ as well as the possibility of discrimination between heavy Dirac and Majorana sterile 
neutrinos via |ALJ_= 2 processes in high luminosity experiments (55| . which is summarized 
in a recent review |33l|. 

Goals: Our goals are the following: 

• a:) Motivated by the possibility of high intensity sources, we assess the signals of 
heavy sterile neutrinos from meson (DAR) (both 7t~;K~) by focusing on searching 
for charged leptons of negative helicity (or positive helicity for their antiparticles in 
7r + ; K + in (DAR)) in a setup akin to the Stern-Gerlach type experiment where opposite 
helicity components are spatially separated by a magnetic field gradient. Meson (DAR) 
produces a monochromatic beam of charged leptons back-to-back with (anti) neutrinos. 
Massive neutrinos yield a negative helicity component for the charged lepton, which 
in a collimated beam may be separated from the (larger) positive helicity component 
by a magnetic field gradient. We study the branching ratio for the negative helicity 
component as a function of the sterile neutrino mass, as a complement to the search 
for monochromatic lines. 

• b:) We assess decoherence effects of sterile-active neutrino oscillations in short baseline 
experiments as a consequence of i) the decay width of the meson, and ii) the stopping 



distance of the charged lepton. As previously found in refs. [56|, |57[ the decay width of 



the meson leads to decoherence of oscillations quantified by the dimensionless ratio 

„ 5m 2 
K 



2ET 



M 



where E is the neutrino energy and Ym is the meson decay width. For example for a 
Pion (DAR) E ~ 30MeV and Y n ~ 2.5 x 10~ 8 eV it follows that K ~ (5m 2 /eV 2 ) and 
there could be considerable suppression of the appearance and disappearance probabil- 
ity in experiments with baseline L ~ 30 — 100 mtsj56|, 57[. Another source of decoher- 



ence is the distance at which the charged lepton is stopped L c : if the charged lepton is 
correlated with the emitted mass eigenstate over a long time scale, the quantum state 
is projected onto an energy eigenstate and oscillations are suppressed j56], |58j]. Both 
effects, meson lifetime and charged lepton stopping scale, are sources of decoherence in 
sterile-active oscillations that are more prominent in short-baseline experiments and 



mass scales 5m 2 ~ eV 2 , as discussed in refs. |56l. |57|. These effects can potentially 
impact the assessment of the sterile neutrino mass, mixing angle and CP-violation 
phases. 

c:) If sterile neutrinos are massive Majorana particles there are neutrino-antineutrino 
oscillations, these are lepton number violating transitions with \AL\ = 2. In short base- 
line oscillation experiments, massive Majorana neutrinos yield two oscillation channels: 
the usual one with AL = and another with \AL\ = 2. While this latter channel is 
suppressed by the ratio m/E, we seek to study these lepton number violating oscilla- 
tions in detail as potential discriminators between Dirac and Majorana neutrinos in 
future high luminosity experiments. Furthermore, neutrino-antineutrino oscillations 
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can distinctly yield information about CP-violating Majorana phases 59[ and one of 
our goals is to assess the impact of the above mentioned decoherence effects on the 
potential measurement of these transitions for new generations of sterile neutrinos. 



Brief summary of results: 

• Massive sterile (anti) neutrinos produced in ir~,K~ decay at rest (DAR) lead to a 
negative helicity (positive if the decay is tt + ,K + ) component of the charged lepton 
produced in the decay. We obtain the branching ratio for this component and argue 
that a Stern-Gerlach filter with a magnetic field with a gradient induces spatially 
separated domains of the different helicity components with a ratio of abundances 
determined by the branching ratio. A combined measurement of the monochromatic 
line for the charged lepton and the ratio of abundances of the spatial domains yield 
simultaneous information on the mass and the absolute value of the mixing matrix 
element. This setup is most sensitive for heavy sterile neutrinos with mass m s in the 
MeV range. We find that for pion (DAR) the electron channel is the most efficient 
for 3 MeV < m s < 135 MeV whereas for K-(DAR) both muon and electron channels 
are similar in the mass range allowed by the kinematics. We obtain an estimate for 
the upper bound on the branching ratio from previous experiments with typical upper 
bounds Br < 10~ 8 — 10~ 6 perhaps accessible in the next generation of high intensity 
experiments. 

• Decoherence suppression in short baseline experiments: the tt, K lifetime V m and the 
stopping or measurement distance scale L c of the charged lepton produced with the 
neutrino introduce energy (or time) uncertainties that lead to decoherence of neu- 
trino oscillations and a suppression of the transition probability. These effects depend 
on the dimensionless quantities 7Z = 5m 2 /2ET M and T M L C . We obtain the gen- 
eral expressions for the appearance and disappearance probabilities that include these 
decoherence effects by implementing a non-perturbative quantum field theoretical gen- 
eralization of Wigner-Weisskopf theory to describe the production and detection pro- 
cesses. We study both Dirac and Majorana neutrinos and show that these processes 
also affect CP-violating transitions. For Majorana neutrinos we study both AL = 
oscillations and | AL = 2| (L is lepton number) transitions from V <H- v oscillations. We 
focus in detail on 3 + 2 and 3 + 1 schemes with new generations of sterile neutrinos and 
obtain the general CP-even and CP-odd expressions for the transition probabilities 
including |AL| = 2 processes with Majorana neutrinos. 

• We analyze in detail the impact of the decoherence effects on current and future 
experimental efforts. The main effects are: i) an overall suppression of the transition 
probability ii) an energy dependent phase shift of the oscillatory terms. We argue 
that if the experimental data is fit to the usual quantum mechanical formula, these 
effects lead to an underestimate of masses, mixing and CP-violating (both Dirac and 
Majorana) angles. 

We find that for vr decay at MiniBooNE/SciBooNE 71 ~ Sm 2 /eV 2 , T M L C ~ 1 and for 
5m 2 > leV 2 , there is a substantial impact resulting in important underestimates of 
the mass of putative new sterile neutrinos as well as the effective mixing angles (for 
the case of 3 + 1 mixing) and CP- violating angles. For example, for a sterile neutrino 
mass m s ~ 3eV we find that the quantum mechanical formulae underestimate the 
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parameters by ~ 15 — 20%. For experiments in which neutrinos are produced by 
the decay at rest (DAR) of ir, K our study suggests that implementing muon (and 
or electron) stopping on distances much smaller than the decay length of the mesons 
would prevent these decoherence effects and the usual quantum mechanical formulae 
would yield the correct fit. We also comment on fundamental differences between 
(current) accelerator and reactor experiments, arguing that the decoherence effects 
are negligible in reactor experiments whereas they affect the transition probabilities 
in accelerator experiments. 



Several appendices provide the technical details. 



II. HEAVY STERILE NEUTRINOS IN RARE n ± ,K ± DECAYS AT REST: 

The possibility of the existence of heavy "sterile" neutrino states had received early 
attention both theoretical lv[6~o| and experimentally (6lM68| . a review of the experimental 
bounds is presented in ref.|69]. In this section we analyze possible observational signatures 
of heavy sterile neutrinos in ir~,K~ — > l~V a decay at rest (DAR) but focus on negative 
helicity charged leptons (or positive helicity for ti + ,K + decay). If the neutrino is massless, 
the charged lepton emerges from tt, K (DAR) with right handed helicity (in the rest frame 
of the meson, which for (DAR) is the laboratory frame). However, if the neutrino is massive, 
a fraction of the charged lepton yield has left handed helicity. If the charged leptons are 
collimated along an axis z and there is a magnetic field that features a gradient along this 
direction, the situation is akin to the Stern-Gerlach experiment: the magnetic field gradient 
leads to a force F oc — V(/2 • B) where /2 is the charged lepton magnetic moment. This 
force spatially separates the charged leptons with spins polarized parallel and antiparallel 
to the magnetic field gradient, just as in a Stern-Gerlach filter. The ratio of the helicity 
population is determined by the branching ratio of the production process into negative 
helicity charged lepton states. Our goal is to obtain this branching ratio, which measures 
the relative intensity of the negative helicity states and could serve as a complement to the 
searches of monochromatic lines. 

While there has been a substantial experimental effort searching for monochro- 



matic lines associated with heavy sterile neutrinos from tt, K decays, we are not aware of 
experimental efforts searching for wrong helicity charged lepton signals in mesons (DAR). 
The bounds obtained from the various experiments[61]-[68] are summarized as exclusion 
regions in ref. [69], which imply mixing angles (rather elements of the active-sterile mix- 
ing matrix) < 10 -6 , making the branching ratios for the se p rocesses very small. However, 
high intensity beams as envisaged in the proposals may provide the experimental 

setting to search for these signals. 

For a 7r or K meson, M, the interaction Hamiltonian for a M — > I v\ decay is given by 



Hi = F M Y, fd 3 x[* la (x,t)^m Va (x, t)J™(x,t)] ; L = i(l- 7 5 ) (II. 1) 

where the label a refers to the charged leptons, J^(x,t) = id^M(x,t) and M is a complex 
(interpolating) field that describes the charged pseudoscalar mesons M = n~,K~. For a 
7r~ meson, we have that F n = V2GFV u dfn and for the K ± meson, we have that Fk = 
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\f2GpV us fx, where are the decay constants. The flavor neutrino fields and the fields 
that create/annihilate neutrino mass eigenstates are related by 

= J2 U °i*»i ■ ( n - 2 ) 

3 

For n generations of Dirac neutrinos the matrix U is n x n, unitary and features (n — 
l)(n — 2)/2 CP-violating Dirac phases. For Majorana neutrinos 

U^U = UD ; D = diag[e iei/2 ,e ie2/2 ,--- ,e^" /2 ] (II.3) 

where [/ is the mixing matrix for Dirac neutrinos and we have allowed an inconsequential 
overall phase. It follows that 

U aj = U aj e^ 12 . (II.4) 
The Majorana CP-violating phases, 9i — 8j, only contribute to v •<-)■ V oscillations and 



|AL| = 2 processes [59j which will be studied in detail in section (|V)) . 

The details of the quantization of the different fields are provided in appendix (jA]). From 
these results, it follows that, after integration over the spatial variables, the relevant Hamil- 
tonian to obtain the production amplitudes is given by (see appendix (TAJ) for notation) 

_ f_m v v v lT & h _y^ = h ' )MM * ~ Mp " } ] k=- - ai5) 
1 ^hbh aj vzem^emem ' p+q ( ] 

where the Fermi quantum fields, if> v ., are expanded as in (jXJ) either for Dirac or Majorana 
fermions. 

We identify the production matrix element M~(p) — » l a {k)V a (q) as 

Ml a (k, q, h, ti) = Uaj M^(k, q, h, ti) (II.6) 

3 

where 

M^(k,q,h,ti) = F M U a!h (k)YLV 3<h iq) Pfl ; k = p + q (II.7) 

is the transition matrix element for meson decay into a charged lepton, a, and an antineutrino 
mass eigenstate, j. For Dirac neutrinos, the spinors Vj t h'(q) m ( 1H.7|) are given by (1A.8I) . 
whereas for Majorana neutrinos 

V j>h ,(q)^Ul h ,(-q) (11.8) 



given by ([Q8jl and the mixing matrix U -^-U given by (1II.3|III.4j) . The separation of helicity 



contributions is frame dependent and the most clear identification of processes that reveals 
a massive neutrino is provided by the decay of the pseudoscalar meson at rest (p = 0) 
so that the laboratory coincides with the rest frame of the meson and helicity states are 
unambiguously recognized. The contributions to the production amplitude from the different 
helicity states in (DAR) are given by 



u a ,4q)i^v j!+ (q) P » = o 
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where (see notation in appendix (|A|) ) 



m n 



E a (q) + q 



VE a (q) + q ! E a (q) = 



; a 



L v 



(11.10) 



Gathering these results, we obtain the helicity contributions to the tt, K decay widths either 
for Dirac or Majorana neutrinos: 



1 w/K-+lv s 



G 



F 

Att 



Ui s \ z \ V ud / us \ 2 fli K q* 



G 



V — " F I TT P IT/ P n * 

L tt/K^Iu s ~ ^ I U ls\ \ v ud/us\ Jjt/K q 



2 r2 



m, 



m. 



EuAq* 



E t (q*) + q* 
Ei{q*) + q* 

E„ s (q*) + q* 



(11.11) 
(11.12) 



where 



2m M 



( m M ~ ( m i + m s) 2 ) (m 2 M - (mi - m s ) 2 ) 



m s < m M - mi 



(11.13) 



and here we refer to the heavy sterile mass eigenstate as s rather than identifying it with a 
fourth or fifth generation. 

In the limit m s — > 0, the usual result for tt, K decay at rest, where the antineutrino and 
the lepton are both right handed polarized, is obtained. We are particularly interested in 
the branching ratio for the process in which both the antineutrino and the charged lepton 
feature left handed helicity, given by (III . 1 2[) . The branching ratio for this process is obtained 
by normalizing to the total meson width and since these are rare processes, we can instead 
normalize to the proxy to the total width 



ytot _ 1 -K/K^jiv 

w/K ~ Br(n/K jiu) 



(11.14) 



where Br(7r/K — > fiu) = 0.999,0.635 is the branching ratio for the purely leptonic decay 
into muons and massless neutrinos for tt, K decay respectively. Specifically, we have 



Br 



r 



~ptot 
1 M 



\Uu 



2 Br(M fj,u) q*m 2 u 



mf,m M 1 3 

M m \ mi 



Ei(q*) + q* 
E Us (q*) + q* 



(11.15) 



Fig. (IT]) show the branching ratios ( III. 151) for tt — > fx, e j7 s and K — > fi, e v s respectively. For 
tt (DAR), the electron channel offers a larger window simply because of the larger amount 
of phase space available whereas the maximum mass available for a heavy sterile neutrino 
in the muon channel is ~ 33.92 MeV. 

A Stern-Gerlach experiment: 

In meson (DAR), the presence of a heavy sterile neutrino is manifest as a monochromatic 
line in the charged lepton spectrum at kinetic energy 



Ti{q*) 



1 

2m M 



(m M - mi) 2 -m 2 ; m s < (m M - m t ) 



(11.16) 
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FIG. 1: Left panel:Br~ 



\Ui s \ 2 , right panel:Br K _ 



\Ui s \ 2 vs. m s for I = fx, e. 



The negative helicity component of the charged lepton in (DAR) provides another mani- 
festation of a massive sterile neutrino which can be exploited in an experiment to complement 
the search of monochromatic peaks in the charged lepton spectrum. The experimental setup 
to exploit the negative helicity component (or positive helicity for the opposite charged me- 
son and charged lepton) should be akin to the original Stern-Gerlach experiment to separate 
spin components. In this case, the relevant quantity is helicity; therefore, consider collimat- 
ing the charged leptons in (DAR) along a z — axis and setting up a magnetic field with a 
gradient along this direction so that the direction of motion of the collimated charged leptons 
coincide with the direction of the gradient of the magnetic field. Under these circumstances, 
there is a magnetic force acting on the charged leptons 



F 7 oc -h- 



dBz 
dz 



(11.17) 



where h is the helicity component; thus, opposite helicity components separate spatially and 
the fraction of negative helicity charged leptons is measured by the branching ratio (III. 151) . 
Therefore, searching for spatially separated domains of charged leptons in combination with 
a monochromatic line, may provide a more robust signature of heavy sterile neutrinos and 
allow extraction of the mixing matrix element \Ui s \: the mass of the sterile neutrino is 
inferred from the peak in the monochromatic spectrum while the ratio of abundances of the 
helicity components is determined by the branching ratio (III. 15ft ; therefore, with the input 
for q* obtained from the peak in the monochromatic line and the measurement of the ratio 
of abundances of helicity states, the branching ratio fill. 15[) yields \Ui s \. 

An estimate for the upper bound on the branching ratios, Br , given by (III. 151) can be 
obtained from the summary of the bounds on the mixing matrix elements, |?7; s | 2 , provided 



in ref . |69| for I 
bound 

and for K (DAR) 



/i: the exclusion region for 7r (DAR) from the fi spectrum yields an upper 
I 2 < 10 5 ; 3 MeV < m s < 33 MeV (11.18) 



U, 



fXS I 



|t^l 2 <io- 6 



10 



-5 



30 MeV < m s < 330 MeV 



(11.19) 
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The experiments 61-68] on which the bounds in ref. [69] are based, search for monochromatic 
peaks in the muon spectrum, both from tt, K (DAR). Assuming that a similar upper bound 
holds for electrons, we find from fig. (pQ) that the upper bound for Br~^^ ei7a 

Br-- flu - s <10- 8 -10~ 7 ; 3 MeV < m s < 33 MeV (11.20) 
Br~- ei7s <10- 8 -KT 6 ; 3 MeV < m s < 135 MeV (11.21) 

The small m s region is obviously suppressed by the m 2 /m 2 factor whereas the region 
near the kinematic edge is suppressed by phase space. For tc decay, the electron channel is 
the most favorable to study the intermediate mass region ~ 3 MeV < m s < 135 MeV with 
typical upper bounds on the branching ratios 10 -8 — 10 -6 . 

For K decay, both fi, e channels yield similar branching ratios with upper bounds in the 
range 

< 9 6 for J 4 MeV < m s < 360 MeV (/i - channel) , , 

nr K^,eu s ~ iu ror < 4MeV < m ^ < 4 14MeV ( e - channel) ' { > 

The "low" mass region of cosmological interest, m s ~ few keV, is much more challenging. 



The experimental results of refs.|6lH68| and the analysis of ref. [69| do not provide reliable 
upper bounds; however, bounds for this mass range emerge from cosmology: a "heavy" 
sterile neutrino can decay into a photon and a light active neutrino, which, for m s ~ keV, 
leads to an X-ray line. Cosmological constraints are summarized in the review articles in 
refs. 42|, |43[ with an upper bound \Ui s \ 2 ~ 10 -10 — 10~ 9 which would make the branching 



ratios exceedingly small, even for the high intensity sources envisaged. 

A firmer assessment of whether the Stern-Gerlach type experiments, combined with 
searches of monochromatic peaks in ir,K (DAR), are feasible in determining the masses 
of "heavy" sterile neutrinos calls for a detailed understanding of backgrounds which is a 
task that is beyond the scope of this article. 



III. OSCILLATIONS IN SHORT-BASELINE EXPERIMENTS 



For short baseline oscillation experiments, the relevant range of neutrino mass differences 
is Sm 2 ~ (eV) 2 . A detailed analysis of oscillation phenomena requires an understanding of 
the production and detection process. In ref. 56[ a quantum field theoretical generalization 
of the Wigner-Weisskopf method [70. 71] was introduced to obtain the correct quantum state 
arising from the decay of the parent particle. A previous treatment of the correlations of 
the decay product within a Wigner-Weisskopf approach to semiclassical wave packets was 
originally studied in ref. 



simple models. In ref. [56 



75 



and the dynamics of propagation were studied in ref.[76j in 
the method was implemented in a simple quantum field theory 
model of charged current interactions and several aspects were found to be much more 
general, such as the decoherence effects associated with the lifetime of the decaying parent 
particle as well as the observation (or stopping) of the charged lepton produced as partner 
of the neutrino in a charged current interaction vertex. 

Many of these aspects were found also in refs. (5?| in a different formulation but without 
explicitly obtaining the quantum mechanical state that describes the decay products. 

Meson decay leads to a correlated state of the charged lepton and the neutrino, a quantum 
entangled stated, H^, 72], the entanglement being a consequence of the kinematics and 
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/s p 

in detail in refs. 56j, |72j, quantum entanglement leads to decoherence in neutrino oscillations 



conservation laws pertinent to the decay[73[. As originally observed in ref.|58[ and analyzed 



which is a result that has been confirmed more recently in [57|,|74| within a different approach. 

In this article, we generalize the quantum field theoretical Wigner-Weisskopf method 
introduced in ref. 56] to describe (pseudoscalar) meson decay via charged current interactions 
in the standard model, including all aspects of the interactions both for Dirac and Majorana 
neutrinos. An alternative formulation is offered in ref. j57[; however, the full quantum field 
theoretical Wigner-Weisskopf method not only illuminates clearly the quantum entanglement 
and correlations between the charged lepton and neutrino states both in momentum and 
helicity, but also allows a systematic study of Dirac and Majorana fermions including the 
dynamics of v <->■ V oscillations and |AL| = 2 processes discussed in detail in section (TVT). 



A. Production from meson decay: 



In appendix |B] (see also ref.[77| for more details), we implement a quantum field theoretical 
version of Wigner-Weisskopf theory and we find the Schroedinger picture quantum state that 
results from pseudoscalar meson decay which is given by 



\Mj(t))) = e-^M'e-r-M ! |M p T(0)) - ]T \ U aj IlJ M^(k, q, h, h>) JF aj [k, f, t] 

q,aj,h,h' 



x e^ EaW+E ^ q))t \lZ(h,k))\vj(h',-^)\ ; k = p + 



Q: 



(III.l) 



where 



n 



V 
aj 



. (III.2) 

[8VE M (p)E a (k)E j (qp 

Although we consider plane wave states, the generalization to wave-packets is straight- 
forward and we comment on the wave-packet approach in section fIVIBj) . The produc- 
tion matrix element Ai^j(k,q,h,h') is given by ( 111.71) . (see eqn. (1B.23I) in appendix (IB1) 
Fm(p) — m M^ m I Em{v) where Y M is the decay width in the rest frame of the meson, and 



J- aj[k, q] t] 



3 -i(B M (p)-s a (fe)-%( g ))t -r M (p)- 



E M (p) - E a (k) - Ej(q) 



(III.3) 



The second term in (IIII.ll) reveals that the emerging charged lepton and neutrino are 
entangled both in momentum and in helicity. 

The factor J- a j[k, q; t] encodes the time dependence of the production process. In order 
to understand the content of this factor, consider the case T M = 0. In this case, 



J- [k, q',t] — e 



-i(E M —E a — Ej) 7t 



2i sin 


(Em — E a — Ej)t 




Em ~ E a — Ej 





t— ¥OC 



-> 2m5(E M - E a - Ej) (III.4) 



namely, in the long time limit, this function describes energy conservation at the production 
vertex. The width of the decaying meson state determines a time (or energy) uncertainty 
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and, either for a narrow width or large time, the function J r a j[k,q;t] is strongly peaked at 
E a + Ej ~ Em which describes approximate energy conservation within the time or width 
uncertainty. 

In a typical experiment, the charged lepton produced by pion (kaon) decay is stopped 
shortly after the end of the pion decay pipe, at which point the correlated quantum state 
after the neutrino state is disentangled by the observation, capture or absorption of the 
charged lepton at t c . 

If the charged lepton l a is observed, or absorbed with momentum k and helicity projection 
hi at time t c , the wave function is projected onto the state (l~ (hi, k)\ and the correct (anti) 
neutrino state that propagates is given by 

\v{q- hi)) = - e -*We J2 U aj Il£ M%(k, q, h u h!) T aj [k, q- t c ) e~ iE ^ \v ,(h! , -q)) , 

3,h' 

(111.5) 

where q = p — k. This neutrino state still carries the label hi as a consequence of the helicity 
entanglement with the measured charged lepton. 

We note that if Ai^j, J^aj, Ej are all independent of the mass of the neutrino, j, these 
factors can be taken out of the sum and the resulting (anti) neutrino state is proportional 
to the familiar Pontecorvo coherent superposition of mass eigenstates. We will analyze 
this approximation below after assessing the total transition amplitude from production to 
detection; however, before doing so, it proves illuminating to understand the normalization 
of the state flTTLo]) . 



UMK) = (V(q ] h i )\V(q ] h i ))= J^K^Ki-Maj&^h^h')] 2 T aj [k,q;t c ] . (III.6) 

j,h' 



In the narrow width limit, the function J-'ajik, Q] t c ] becomes oc 5(E M (p)—E a (k)—Ej(q)) 
and the proportionality constant can be obtained by integrating this function in the variable 
£ = Em(p) — E a (k) — Ej(q), from which we find 



J~ otj\k, q, t c 



2tt 



Tm(p) 



I _ e -^M{p)tc 



S(E M (p)-E a (k)-E 3 (q)). 



(III.7) 



Therefore 



N v (q] K) 



Tm(p) 



\ — e -Thl{p)tc 

' J2\ U ^\ 2 K j Ml J (k,q,h l ,h')\ 2 2nS(E M (p) - E a (k) - E 3 (q)) . 
3,h' 

(111.8) 

In appendix ([C]), we obtain the relation between the normalization ( 1III.8j) . the partial 
and total decay width of the meson and the number density of charged leptons produced 
by meson decay during a time t c . While ref. 13] discusses the normalization of the neutrino 
state 1 , to the best of our knowledge, the relation of the neutrino normalization to the number 
density of charged leptons produced has not been recognized previously. 



1 See section (8.1.1), pages 285,286 in ref. [13 1 . 
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B. Detection via a charged current vertex: 



In what follows we assume the neutrino to be described by a Dirac fermion, extending 
the discussion to Majorana fermions in section (jVj) . 

Consider the case in which the (anti) neutrino is detected via a charged current event 
V N — > 1 1 N' at a detector situated at a baseline L (fig.©). 




N 

FIG. 2: Production via M~ — > l~V detection via a charged current vertex V N — > it N' at a baseline 
L with N, N' nucleons or nuclear targets. The charged lepton l~ produced with the antineutrino 
is observed, absorbed or decays at a time t c . 



The Schroedinger picture quantum states that describe the initial and final states are 



\i) = \V;N) = \v)®\N) v ; |/) = |Z+;JV , ) = |Z+)<g»|JV , > 



(III.9) 



where \v) is given by flIII.5j) . the state \N)x> describes a nucleon or nuclear target localized 
at the detector and the outgoing charged lepton is measured with helicity hf. The transition 
amplitude in the Schroedinger picture is given by 



(f\e 



-iH(t n -to) 



i) 



-ie 



iEptD [ ° e iEFt ' (f\H ie - lHot ' e iHotc \i) dt' (111.10) 

Jtr 



where E F = + E N > is the total energy of the final state, and H , Hi, H are the unper- 
turbed, interaction and total Hamiltonians respectively. To obtain this expression we have 
used e-^^o-tc) _ e -iH t D u(t D ^ e iH o*c anc [ jj(t D ,t c ) is the usual time evolution operator 
in the interaction picture. 

Up to an irrelevant overall phase we find 



T^f = KP ai M v a] J~ a j Gpj (1+ ; N'mv^ ; N) 

j,h' 



(III.ll) 



where 11^ is given by eqn. 2j) . we have suppressed the indices to avoid cluttering the 
notation and introduced 



G 



Pi 



e 2 * 



2 sin 

Ei)(t D +t c ) 


\{E F - 


E N — Ej)(t D — t c ) 




Ep 


— E N — Ej 





(111.12) 
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The relevant interaction Hamiltonian is given by 



Hi = U; 3 V2G F J ^.(f) 7 M L ^ h {x) j( N ' N '\x)d 3 x + h.c. (111.13) 

where J^ N ' N ^ (x) is the hadron current with matrix element 2 

jN,N'(p\ 

(N\J^ N '\x)\N') ee ]T feM eiP ' S - ( IIL14 ) 



(1+ ; N'lH^V, ; N) = U* pj nj, MZ (HI.15) 



leading to the matrix element 
where 

n?,- = r (HI.16) 

^ = V^G F V i , ft /(g)7'']LV i 8, h/ (^)^(P) ; q = k-p = P + k'. (111.17) 
Therefore, the total transition amplitude from production to detection is given by 

7^/ = £ u aj a<J ^ g a - nj ; (in. is) 

where we have suppressed all the arguments to simplify notation. The factors T a j and Gpj 
encode the time dependence of the production, measurement of the charged lepton produced 
with the (anti) neutrino and final detection processes and the energy uncertainty from the 
finite lifetime of the parent meson. As noted above (see eqn jIII.4l ). JF a j describes nearly 
energy conservation in the long time narrow width limit but includes the energy uncertainty 
from the width of the decaying state. Similarly 

Gpj'^F 2-k8(E f - E N - Ej) (111.19) 

describes energy conservation at the detection vertex in the long time limit. 

The phases in these factors encode the information of interference effects between the 
different mass eigenstates. 

In order to isolate the contribution of these factors there are several approximations that 
are dictated by the experimental aspects: 

Approximations: 

1. For neutrino masses consistent with baselines of few hundred meters, namely rrij ~ eV 
and typical neutrino energy from meson decay > 30 MeV, neutrinos are ultrarelativis- 
tic so we can approximate Ej(q) = E(q) + m 2 J2E(q) with E(q) = q. 



2 This matrix element may be written in terms of vector and axial vector form factors, but such expansion 
is not necessary in our analysis. 
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2. We neglect the neutrino masses in the factors Ej(q) in the denominators in IT^. , IT?- 
feqns JIlX2TITLT6|) . 



3. We also neglect the mass dependence of neutrino spinors V, which depend upon the 
mass through the factor Sj(q) = rrij/(Ej(q) + q) (see (IA.8|IA.9j) ). Neglecting the neu- 
trino masses in the spinors leads to the production and detection matrix elements 
J\A V , M v to be independent of the neutrino masses, therefore independent of the 
label j. 

4. Neglecting the neutrino mass, the negative chirality (anti) neutrino only features a 
positive helicity component, therefore only h! — + remains in the sum. This is, 
obviously, a consequence of £j < 1. 

Under these approximations and the unitarity of the mixing matrix U, the normalization 
N v (q) (1III.6j) of the neutrino state \V(q)) becomes 

|2 



I — g- r M(p)*c 



2tt5(£ v ) ; £ V = E M -E 



E. 



(111.20) 



\U v M vl 
Tm(p) 

The time dependent factors G@j feature phases whose interference leads to the oscil- 
lations in the transition probabilities, therefore the terms m 2 -/2E{q) must be kept in these 
phases. 

Under these approximations, the factors Tl v M. v and U V A4 V can be taken out of the sum 
and the final result for the transition amplitude factorizes into production, propagation with 
oscillations, and detection contributions: 



U P M V 



03 



(111.21) 



Production Propagation-Oscillations Detection 



The transition probability is given by 



Pi 



n? M 



V 



(111.22) 



where T a j = T a j\k, q, t c ] is given by (1III.3j) evaluated at t — t c and Gpj is given by (1III.12j) . 



It proves convenient to introduce: 

£ v = E M {p) - E a {k) - E{q) 



& 



v 



Ep — E 



N 



E(q). 



(111.23) 



For rrij <C E(q) and narrow width Tm <C Em, the products T a jT^i are sharply peaked at 
£-p, becoming nearly energy conserving delta functions in the long time and small width 
limit (see f THI.4j) ). Similarly, GpjG*^ is sharply peaked at E-p- Each term J 7 , G describe 
approximate energy conservation at the production and detection vertices respectively. In 
order to extract the coefficients of the energy conserving 5(£p), 8(£t>), we integrate the 
respective products with a smooth initial and final density of states that are insensitive to 
T M and Aj (for details see ref. [56|). We find 



J Q-3 J ai 



2tt 



f(p)* c 



Tm(p) 



1 + iKi 



8{£ v ) 



(111.24) 
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were we have introduced 



A -. 6m i ■ v A v 

ij o 771/ \ i ij 



5m% E M (p) 



2E(q) 



similarly 



T M (p) 2T M M M E(q) 

\ _ g*A<j(txj— -c) 

G Pj G% = 2vie i ^ 



bm ii = m i — 



A,- 



(111.25) 



(111.26) 



d 
dtn 



As usual, one is interested in obtaining the transition rate; therefore, we focus on 
7I_>/| 2 for which we need 



dt 



D 



(111.27) 



Separating the diagonal % = j from off-diagonal terms in the sums in 22|) . and using the 
result ( 1IH.20|) for the normalization of the (anti) neutrino state, we find the transition rate 



dt 



\% 



D 



v, 



dr 



(2n) e V 2 d 3 k'd 3 P 



where 



2tt 6(£- 



(111.28) 



(111.29) 



(27v) e V 2 d 3 k'd 3 P 

is the double differential detection rate for VN — > ItN' for an incoming massless neutrino, 
and V a ^/3 is the flavor transition probability 

V a ^ = U^TPJUpila (111.30) 



0,1 



where 7jj are the interference terms 



1 _ e -*Ayt- e -r M (p)t- 



1 _ e -r M (p)* c 
Unitarity of the U matrix allows to write 

7V>/3 = -W» - 2^Re[ 



1 — ilZij 



i = i* 



(111.31) 



2 

j>i 

-2 

j>i 



UajUpjUaJJpi 



Re 



1-1. 



Im 



1 y 



(111.32) 



In the above expressions we have implicitly assumed Dirac neutrinos, the case of Majorana 
neutrinos is obtained by the replacement (see eqns. (jII.3)II.4p ) U — > U ; U a j = U a j e 1 ^^ 2 Va 
from which it is obvious that the CP-violating Majorana phases do not play any role in 
V a — > Vp oscillations. 

The possibility of CP violation in the neutrino sector from Dirac phases is encoded in 
the imaginary part in (lIII.32j) since for the transition probabilities for u a — > z/g it follows 
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that Utf — > U* { . Therefore decoherence effects in the imaginary part of Jy lead to possible 
suppression of CP-violating contributions in the transition probabilities. 

The transition rate (1111.280 . along with ( 1111.320 . are some of the important results of 



this article; the factorized form of ( 1111.280 is a consequence of the approximations described 
above. The origin of the prefactor J\f v is clear, it is the normalization of the neutrino state 
that emerges from disentangling the charged lepton in the production process, since this is 
the correct neutrino state that propagates to the detector and triggers the charged current 
reaction that yields the measured charged lepton in the final state. The interference terms 
( 1IIL31I) encode the decoherence effects arising from the finite lifetime of the source and the 
energy uncertainty associated with the time scale in which the charged lepton produced in 
a correlated quantum state with the neutrino is observed (or captured). This decoherence 
can be understood clearly in two limits: 

• When Aij <C T M it follows that IZij — > and ly is the usual interference term. 
In this limit the energy uncertainty associated with the lifetime of the source does 
not allow to separate the mass eigenstates and the coherence of the superposition of 
mass eigenstates is maintained. However, in the opposite limit, A y - ^> T M , the factor 
TZij ^> 1 and the interference term is suppressed. In this limit, the lifetime of the 
source is long, the corresponding energy uncertainty is small and the mass eigenstates 
are separated in the time evolution and coherence between them in the superposition 
is suppressed. 

• In the limit T^j — > it follows that 

j e <A„(to-W2) sm[A tJ t /2] (m } 

There are two effects in this expression: 1) a shortening of the baseline by the distance 
travelled by the charged lepton produced with the (anti) neutrino and 2) a suppression 
factor associated with the time uncertainty: if Ajji c > 1, then the interference term is 
suppressed, this is because if the charged lepton produced with the (anti) neutrino is 
entangled all throughout the evolution at long time t c 3> l/Ay the energy uncertainty 
becomes much smaller than the difference in energy between mass eigenstates and 
these are projected out by energy conservation which leads to their decoherence in 
the superposition. This is another manifestation of energy conservation as encoded in 
Fermi's Golden rule. In terms of the oscillation length, L°? c , defined as 

5mf, 2n . 
A w = — ^ = 111.34 

the suppression factor 2 sin[Ayt c /2]/Ayt c < 1 when the stopping length scale L c = 
t c ~ Lij. The case r — > is relevant for reactor experiments. See the discussion in 
section(|VL£}. 

The suppression factor associated with the lifetime is relevant in the case of possible new 
generation of (sterile) neutrinos with masses in the eV range when produced in the decay 
of pions or kaons. 

For pion decay at rest, the typical energy of a (nearly massless) neutrino is E* ~ 30 MeV, 
the pion width at rest = 2.5 x 10 _8 eV and for one generation of sterile neutrino with 
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m 4 ^> mi ]2i3 we find 

Kc^c 2 -(^±)\ (111.35) 

therefore, for m 4 > 1 eV, the suppression factor can be substantial and the transition proba- 
bility is suppressed. For the decay of a pion in flight with a large Lorentz 7 factor, the result 
only changes by a factor 2 as can be seen as follows: consider a neutrino that is emitted 
collinear with the direction of the pion in the laboratory frame (say along the z — axis), its 
energy in the laboratory frame is 

E = -fE*(l + V n ) (111.36) 

where V n is the pion's velocity, for 7 ^> 1 it follows that E ~ 2jE*. The width of the pion 
in the laboratory frame is r^/7; therefore, for neutrinos produced by pion decay in flight 
with a large Lorentz factor 

In conclusion, for new generations of (sterile) neutrinos with masses in the eV range, 
experiments in which oscillations are probed with neutrinos from pion decay feature the 
suppression factors associated with the pion width. For Kaons, the situation improves 
because in this case 

~ 5 x Kr 8 eV ; E* ~ 235.5 MeV 

and for Kaon (DAR) 



K ~ 25 KeVJ ' 



thus 1Z < 1 for m 4 ~ few eV. 



IV. 3 + 2 AND 3 + 1 CASES IN THE "SHORT-BASELINE APPROXIMATION": 



In the "short-baseline" approximation, we assume that there are sterile neutrinos j = 
4, 5--- with m 4 ,m 5 --- 3> m 1 ,m 2 ,m 3 so that 5m 2 L/E ~ (9(1) for L ~ 10 — 1000 mts 
corresponding to short baseline experiments. 

We begin by considering the 3 + 2 scenario from which we will extract the case 3+1. 

3+2 case: In this case, m 5 ~ m 4 ^> mi,m 2 ,m 3 so that 

— 1 ,i,j = 1,2,3 ; IiA = hA ; I i5 = h 5 , i = 1, 2, 3 . (IV.l) 

Unitarity of the U matrix entails 

3 

KP* = 5 *e - KiUi* - u^U/k . (iv.2) 

i=l 
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Separating the terms with j = 4, 5 in (1III.32I) . we find for a ^ ft (appearance) 



Mu a4 \\u p4 \ 



\U a i\\Up4\ 



U a5 \\U^ 5 \ COS ( 



'51 



'41 



- 4|C/ a4 ||C/ i 34||L' Qi5 ||C/ i g 5 | cos054 - Re[l - J 54 ] 



+ MU a5 \\U p5 \ 



\Ua5\Wp5\ 



U a4 \\Up±\ COS i 



'54 



Re[l 



'51 



+ 2 



U a4 \\U 04 \\U ari \\U^\ sin, 



'54 



Im 



I 4 i - I51 + I 



L 54 



where following 



we have defined 
054 = Arg 



U a5 Up 5 U a4 Up4 



for a = e , ft = // 



(IV.3) 



(IV.4) 



and used Im[iy] = — Imfi^]. We note that interchanging a •<->• ft (e +»• //) is equivalent to the 
exchange 4 H 5, namely 54 — >■ —^54 = ^45 which leaves the result flIV.31) invariant since 
Ref/jj] is even and Imlji odd respectively under i +»• j. If 54 7^ 0, there is CP- violation in 
the neutrino sector because <p^ 4 — > —<p^ 4 for v — > v oscillations since this implies that the 
elements of the mixing matrix U ai — > U*^ 

The 3 + 2 case effectively describes mixing between three species; consequently, it features 
only one effective CP-violating angle. 

For a = ft (disappearance), we find 



1-M \U, 



ol 



l-\U„ A \ 2 -\U, 



a5 



-Re[l-J 41 ] 



+ |f^«4| 2 |^a5| 



Re[l - I 



54J 



+ \u, 



a5 I 



1 - \u a4 \ 2 - \u, 



a5 I 



Re[l - J 51 ] 



(IV.5) 



which does not feature a contribution from the CP-violating angle. 

3+1 case: This case is obtained from the 3 + 2 case above by setting U a § = Va, leading 
to the appearance probability (a 7^ ft) , 



1 



— 4|[/ Q 4| \Uj34 

and the disappearance (survival) probability (a = ft) 



2 ^Re[l - / 



4lJ , 



7W = l-4|C/ ( 



a?4 



I" If/, 



a 1 



Re[l 



'41 



(IV.6) 



(IV.7) 



The 3 + 1 case effectively describes mixing between two generations and, consequently, 
does not feature any CP- violating contribution. For this case, it is often convenient [l3j to 
introduce the effective mixing angles 



sin 2 2# Q/3 = 4\U a4 \ 2 \Ufs 4 \ 2 , a^ft 
sm 2 28 aa ee A\U a4 \ 2 1 - \U a4 \ 2 , a = ft . 



(IV.8) 
(IV.9) 
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V. MAJORANA STERILE NEUTRINOS AND |AL| = 2 v O V OSCILLATIONS: 

In the previous section we have assumed that sterile neutrinos are of the Dirac variety; 
however, if neutrinos are Majorana fermions, new processes, such as neutrino-less double 



beta decay (see[78[ for recent reviews) and v -h- v oscillations, are available. As discussed in 
ref. [59|], v V oscillations have the potential to reveal CP- violating Majorana phases and, 
to make clear the Majorana nature of the mixing matrix, we write it as U following eqn. 

(HH). 



These processes can be understood by considering the full interaction Hamiltonian in- 
cluding the the hermitian conjugate of the one displayed in flIII.13j) . namely 

+ U Pj V2G F J ^ h {x)r^^{x) J^ N ' N '\x)£x . (V.l) 

The first line yields the AL = V <H- V oscillations just as for the Dirac case discussed in the 
previous section. The second line contributes to the detection process only for Majorana 
neutrinos and yields the |AL| = 2 contribution, as can be simply understood from the 
following argument pertaining to Majorana fermions: the production Hamiltonian flll.lj) 



is determined by charge conservation: a 7r decays into a negatively charged lepton l a , 
thus requiring the \f^ Q in ( IH.ip . the creates a neutrino (same as an antineutrino for 

Majoranas) with an operator bt h that multiplies a charge conjugate spinor U^(k) (see the 



expansion ( 1A.17|) in appendix [A]) . Using the first line in fIV.lj) . the neutrino is destroyed at 
the detection vertex using a b^ h of the which also multiplies the spinor U£{k) along with 
the creation of positively charged lepton It. Therefore, this AL = contribution is the same 

as that for a Dirac neutrino and features the product U a jU^, which is insensitive to the 
Majorana phase. However, the neutrino in the intermediate state can also be annihilated 
by using a from ty v in the second line, which now multiplies the spinor Uh{k), along 
with the creation of a negatively charged lepton lg from \1/ ' J g. This contribution features the 
product U a jUpj = U a jU/3j e %di and manifestly displays the Majorana phase. This process is 
depicted in fig.(E])- 

The CP-conjugate process 7r + — > /+ v — )■ uN — > N' 1^ with | AL\ = 2 features the product 
U*jUpj showing that the Majorana phase is also CP- violating. It is convenient to introduce 

3f = (W) (V.2) 

we now find for the transition matrix element 71-*./ from the initial (i = tt~ N) to the final 
(/ = laljN') state 



V2F M G F E VaiV M (Kip) l (a ■ J M ) (U^{q)) £ JF aj H^ 



h j 



x (ul h „) L (a ■ J^) (u h M)G M Til (V.3) 

where again we suppressed arguments to simplify notation. The sum over helicity states h 
can be carried out straightforwardly using the results of appendix ((A]), we find (no sum over 
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N 

FIG. 3: |AL| = 2 process from Majorana neutrinos. The charged lepton l~ produced with the 
neutrino is observed, absorbed or decays at a time t c and another charged lepton 17 is detected. 



a, (3) 

where we have introduced 
T a f = V2F M G F 2E(q)UlU 



nap 



na/3 



Y, u « u > 



Pi 



2E ] {q) 



(V.4) 



v 



a-J M )v a (q) 



U h,h' 



a 



; a, b 



(V.5) 

where the Weyl spinors v±(q) are the helicity eigenstates f ]A.12j) . Here, E(q) = q for massless 
neutrinos and T°£ do not depend on the mass eigenstate label j. In arriving at expressions 
f lVT4TV3|l we have written I t? rfc = Iff n£ 2E(q)/2E j (q) where the Iff now correspond 
to the phase space factors ( 1111.21111.161) for massless neutrinos, namely Ej(q) — > E(q) = q. 

The amplitudes T®f have a simple interpretation: T"^ is the amplitude for the combined 
process ir~ — >■ /~z/ , z/iV — >■ iJiV' and T^! for the process 7r~ — > l~v , z/iV —> IgN' where u, V 
are massless left handed neutrinos (and right handed antineutrinos) , corresponding to the 
v •<->■ v mixing that violates lepton number by two units. These two amplitudes contribute 
coherently to the process tt~N — > l~l^N' and are added (with their respective signs) in the 
total amplitude. 

The expression ( IV. 4 p is generally valid for arbitrary masses of Majorana sterile neutrinos 
and, with a simple modification of the final state, also describes the AL = 2 processes 



studied in ref. 55 



Proceeding as in the AL = case, we finally find for the transition rate 



d 



dt 



AL\=2 



D 



no sum over a, (3 



(V.6) 



where 



P^p 2 — UgjUpjUgjUm 



3,1 



mjirii 



(V.7) 
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is the v ■<->■ v transition probability with |AL| = 2 and 



T 



T 



a/3 



T 



a/3 



2n5(£ v ) W% 2tt5(£- 



v 



(V.8) 



r M (p) L 

encodes the transition matrix elements for production and detection. We note that unlike 
the AL = case, here there is no factorization of production and detection, this is a 
consequence of the fact that v -h- V oscillation implies helicity change (and a mass insertion) 
and both helicity changing contributions contribute coherently to the total amplitude as 
explained above. A similar observation was pointed out in ref.j79|. We are not concerned 
here with T a/3 but with the transition probability Pjfg L ^~ 2 , which can be written as 



D \AL\=2 



,2 m ) 



j>i 



AEjEi 
mjirii 



+ 2 £ Mu^MKm 4^ Mia] . (v.9) 



j>i 



Just as in the AL = case, the main difference with the usual quantum mechanical case 
is the replacement 



1 _ e ~iAijL c e -T M (p)Lc 



1+n 



(V.10) 



where Ay; 7^.^- are given by eqn. dIII.25|) where the extra factors describe decoherence effects 
associated with the lifetime of the decaying meson and the measurement of the charged 
lepton partner of the produced neutrino. 

To be sure if the absolute mass scale of the new generation of sterile neutrinos is ~ eV 
then the factor ~ m 2 /E 2 < 10 -14 makes the |AL| = 2 contribution all but unobservable 
with the current (and foreseeable) facilities for short-baseline experiments with m ~ eV. 
However, oscillation experiments measure the squared mass differences; therefore, in absence 
of a determination of the absolute scale of masses, there remains the possibility that new 
generation of sterile neutrinos may be heavy but nearly degenerate so that the difference 
in squared masses is small and lead to interference and oscillations on the length scales of 
short baseline experiments and Pj^p~ 2 is not negligible . 

3+2 and 3+1 schemes: Under the assumption that 777.4,777.5 ^> rrii,i = 1,2,3, the con- 



tribution from active-like mass eigenstates is clearly sub leading for the |AL| 
therefore, keeping only the two largest mass eigenstates 



2 transitions; 



P 



\AL\=2 



u a5 \ 2 \u, 



B5\ 



,2 m l 



AEi 



u aA v \u, 



34 1 



12 m l 



4Ej 



+ 



2|C/a5||^35||^a4||f//34| 003(554 + 6*. 



5 1 1 



+ 2|C/ a5 ||L^ 5 ||C/ Q ,4||^/34| sin(<5 5 4 + 54 ) 



m 5 m 4 
AE^El 



Re[/ 54 ] + 
Im[J 54 ] + 



(v.n) 



4E 5 E 4 

where the dots stand for the contributions from % — 1, 2, 3, U is the Dirac mixing matrix 
(JID| and 
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is a Dirac CP-violating phase different from the 054 that enter in the AL = case (1IV.4j) 
and 6*54 = 9 5 — # 4 with 9j the Majorana CP-violating phases (111.31) . 

The 3 + 1 scheme is obtained simply by setting U a 5 = OVa in which case there are no 
oscillations to leading order in m/E. 



VI. ANALYSIS OF DECOHERENCE EFFECTS IN ACCELERATOR EXPERI- 
MENTS: 



The decoherence effects associated with the lifetime of the source and the measurement 
(or capture) length scale of the charged lepton emitted with the (anti) neutrino are encoded 
in the quantities Re[Jjj] , Im[i^] given by eqns. ( lVI.3|VI.4p the latter one determines 
the suppression of the CP violating contributions from these decoherence effects. In this 
section we compare these terms to the familiar ones obtained from the quantum mechanical 
description of neutrino oscillations (IVI.6j) as a function of the neutrino energy for fixed 
baselines. 

Introducing 



Aji(E) 
and replacing as usual 
we find 
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we note that 



(VI.5) 



this relation highlights that there are only two combination of parameters that determine 
the corrections, namely TZji and T M (p)L c ; furthermore, T M (p)L c = L c /l M (p) where Im{p) is 
the decay length of the meson in the laboratory frame. 

In absence of the decoherence contributions, the usual expressions emerge, namely 



Re [J. 
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2E 



sm 



2E 



(VI.6) 
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with 



5m% L 1 5m% 

— 3 — = 2.54: 
2E 



L\ ( GeV^ 



(VI.7) 



eV 2 I V km J \ E 

Whereas the length scale L c is determined by the particular experimental setting and is 
therefore a parameter, the width of the parent particle is a function of the neutrino energy 
through the Lorentz factor as follows. 

In the rest frame of the decaying meson, its width is F m and the antineutrino (neutrino) 

is emitted isotropically with an energy E* = J q* 2 + m 2 with q* given by flII.13j) ; in the 

laboratory frame, where the meson is moving with velocity Vm, the width is Vm/i and the 
energy of an anti (neutrino) collinear with the meson is blue shifted to 



E = 1 E*(l + V M ) 
where we have neglected the mass of the neutrino. Therefore 

E 2 + E* 2 
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In the analysis below, we focus on neutrinos from Pion decay and the analysis for Kaon 
decay is similar. Using the Pion decay width, = 2.5 x 10~ 8 eV, as a benchmark, we 
obtain 



„,„, ( 30MeV \ E- 2 

An illuminating interpretation of the results flVI.3|VI.4p emerges by defining 3 
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in terms of which we find 
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Note that sign^j) = sign^m^). 
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While the general case must be studied numerically, the limit TmL c 3> 1 provides a most 
clear assessment: as compared to the usual quantum mechanical expression (1VI.6j) . the 



decoherence factors result in i) a suppression of the transition probabilities ~ 1/ yl + 7^ 2 j 

and ii) an overall energy dependent phase shift 6ji(E). 

For example, for a sterile neutrino mass m s > 1 eV 3> ^1,2,3 from ir decay, it follows that 
1 < 1Z thus from (1VI.ll|VI.12j) tt/4 < 8ji < 7r/2. Trying to fit the mass (and mixing angles) 



by using the usual expression (1VI.6j) would imply an effective <5^ 2 // = — 2E 8(E) / L. 
For example, for accelerators experiments with E ~ GeV , L ~ 1 Km, such a fit would lead 
to IE 8(E) I L ~ eV 2 and a /arge underestimate of the sterile neutrino mass and the mixing 
and CP-violating angle. 

A similar interpretation holds for the imaginary part ( lVI.14p . which is associated with 
CP-violating amplitudes, the suppression factor would result in an underestimate of CP- 
violation if the usual quantum mechanical expression flVI.6j) is used in fitting experimental 
data. 

This simpler case illustrates that for short baseline accelerator experiments in which neu- 
trinos are produced from the decay of pions and are designed to reveal oscillations of new 
generations of sterile neutrinos with masses in the eV range, the decoherence aspects as- 
sociated with the pion lifetime and the stopping length scale of the muon comparable to 
the decay length of the pion may lead to substantial corrections to the quantum mechan- 
ical oscillation probabilities. A more reliable assessment is obtained numerically below for 
different experimental situations. 

MiniBooNE/SciBooNE: For MiniBooNE/SciBooNE, antineutrinos are produced pri- 
marily from 7r~ — > fi~V^, Pions decay in a decay tunnel ~ 50mts long and muons are 
stopped in the "dirt" at a typical distance ~ 4mts beyond the decay tunnel 4 , therefore in 
this situation T n (p)L c ~ 1. The SciBooNE detector is at a distance L = 100 mts from the 
production region, in between the end of the decay tunnel and MiniBooNE, whose detector 
is at a baseline L = 540 mts, and the neutrino energy range (for both) is 0.3 < E < 1.6 GeV. 
Figs. (jUEJ) show the comparison between the CP-even and odd parts with (modified) and 
without (QM) the decoherence corrections for MiniBooNE for m s = l,2,3eV respectively 
and figs. (EJUS) show the same comparison for SciBooNE parameters with L = 100m and 
same energy range and values of m s . 

These figures confirm the interpretation of the decoherence modifications in terms of an 
overall suppression of the amplitude and a phase-shift that leads to an offset in the position 
of the peaks with respect to the quantum mechanical result. Since the mixing angle is 
extracted from the maximum amplitude of the probability and the mass from the position 
of the peaks, a fit with the quantum mechanical formula would underestimate both the 
mixing angle and the mass as analyzed above. A similar conclusion applies to the CP- 
violating angle. For MiniBooNE, the suppression and off-set are small when 5m 2 < leV 2 , 
resulting in an underestimate of about 3 — 5% in amplitude and mass as shown in fig. (J3|) 
but is larger at SciBooNE as shown in fig. ([7]), but for 5m 2 = m 2 ~ 3eV 2 , fig. ([6]) for 
MiniBooNE reveals ~ 15% suppression in the amplitude with a similar underestimate in the 
mass (off-set). 



Pions and Kaons (DAR): Recent proposals [52|, |53j for high intensity sources to study 



sterile-active oscillations with pions and kaons (DAR) motivate a study of the decoherence 



D.B. is indebted to William C. Louis III for correspondence clarifying these aspects. 
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FIG. 4: CP-even/odd parts of transition probability for MiniBooNE parameters: L = 
540m, T n L c ~ 1 for m s = leV. Solid line (modified) Re[l — I s i]/2 dashed line (Qm) is the 
quantum mechanical result sin 2 [m^/AE]. 




FIG. 5: Same as Fig.© for MiniBooNE for m s = 2eV. 



effects in these experiments. For (DAR) the energy is fixed at E = E* and presumably the 
baseline L is also fixed, we take L = 30 m as a middle-range indicative value for the purpose 
of our analysis, other values can be explored numerically. What is less clear is the value of 
the product T M L C which will ultimately depend on the experimental design. We study the 
possibible ranges TmL c < 1,~ 1,> 1 respectively as a function of m s . For tt — K (DAR) 
it follows that E* = 29.8 MeV ; E* K = 235.5 MeV respectively for which we find the ratio 
f IVLTTj) to be 

K^o = f($) ; «^) = s(^) < VU5 > 

The comparison between the modified results ( 1VI.3IIVI.4I) and the usual quantum mechan- 
ical results ( 1VI.6I) are displayed in figs. fllOH15p for tt, K (DAR) for a baseline representative 
L = 30 m. 
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FIG. 6: Same as Fig.® for MiniBooNE for m s = 3eV. 
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FIG. 7: CP-even/odd parts of transition probability for SciBooNE parameters: L = 100m, T n L c ~ 
1 for m s = leV. Solid line (modified) Re[l — / s i]/2 dashed line (Qm) is the quantum mechanical 
result sin 2 [m^/4E]. 



It is clear from this analysis, both for tc, K (DAR), that decoherence effects are very 
small and the modified result is indistinguishable from the usual quantum mechanical results 
( 1VI.6I) whenever TmL c <C 1 but become large for TmL c > 1. The decay length for tt, K are 



1-x = 7.8 m, Ik = 3.7 m respectively; therefore, in order for the usual quantum mechanical 
results ( 1VI.6P to describe a correct fit to the experimental data, the design must ensure that 
charged leptons (mainly //) be stopped at distances L c 1-k,Ik respectively, namely a few 
cm beyond the stopping target of the mesons. 

Long baseline experiments: For long baseline experiments, the decoherence terms do 
not contribute. This is because these experiments study oscillations with 5m 2 ~ 10 _3 eV 2 
and E ~ few GeV for which L ~ 300 — 1000 km, an example of such experiment is Minos 
in which pions produce neutrinos as in MiniBooNE/SciBoone. In these experiments, 1Z < 
10~ 3 ; TmL c < 1 so that AijL c 1; therefore, decoherence effects are all but negligible 



26 




FIG. 8: Same as Fig. © for SciBooNE for m s = 2eV. 




FIG. 9: Same as Fig. £T} for SciBooNE for m s = 3eV. 



generally for long baseline experiments. 



A. Reactors vs. accelerator experiments 



The suppression of the transition probabilities through the decoherence effects depend 
both on the lifetime of the parent particle and the stopping distance of the charged lepton 
which is produced along with the (anti) neutrino via the charged current interactions. This 
establishes a, fundamental difference between accelerator and reactor oscillation experiments: 
whereas in accelerator experiments neutrinos are produced via the decay of short lived 
mesons with typical lifetimes ~ 10 _8 secs and widths ~ 10 -8 eV, in reactors the (anti) neutri 



nos are produced via the (3 decay of long-lived unstable nuclei 235 U , 238 U , 239 Pu , 241 Pu [2 f 1.122 



with typical lifetimes in the range between hundreds and thousands of years. Furthermore, 
in current short baseline accelerator experiments such as MiniBooNE/SciBooNE, pions de- 
cay in a decay pipe and muons are stopped at short distances beyond the decay pipe so 
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FIG. 10: 7r (DAR) CP-even and CP-odd contributions for T n L c = 0.01 vs m s . The solid line 
(modified) shows Re[l — I s i] where Re [/si] is given by (|VI.3|VI.4p . the dashed line is the quantum 
mechanical result (|VI.6P for Sm^ = m^. 



71 DAR, TU-1, L-30m 




TtDAR, rL c -l,L-30m 




0.0 0.5 



FIG. 11: Same as fig. {TD]) with T V L C = 1. 



that T^Lc ~ 1; in reactor experiments, muons are stopped in the reactor core on distance 
scales so that TL C , AjiL c 1. Our study above clearly shows that under this circum- 
stance the modifications from decoherence are negligible and the transition probabilities 
are indistinguishable from the quantum mechanical result. Therefore, we conclude that the 
quantum mechanical fit to the oscillation probabilities in reactor experiments is always jus- 
tified, whereas in accelerator experiments, decoherence effects both from the lifetime of the 
parent meson and the stopping length scale of the charged lepton partner are substantial for 
m s > 1 eV and a fit to the usual quantum mechanical transition probabilities both for CP- 
even/odd contributions may substantially underestimate masses, mixing and CP-violating 
angles. 



28 



II DAR, TLc-100, L-30m 71 DAR, rL -100, L-30m 




0.0 0.5 1.0 15 2.0 2.5 3.0 3.5 4.0 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 



m,(eV) m,(eV) 




FIG. 13: K (DAR) CP-even and CP-odd contributions for TkL c = 0.01 vs m s . The solid line 
(modified) shows Re[l — I s i] where Re[/ S i] is given by (|VI.3|VI,4p . the dashed line is the quantum 
mechanical result (|VI.6P for 5m? sl = ml. 



B. Wave packets: 



Our study is restricted to plane waves to exhibit the main results and conclusions in 
the clearest possible setting. As has been argued in the literature 8^ 87 ], wave packet lo- 
calization may be an important ingredient in the description of neutrino oscillations. The 
localization length both of the production and detection regions define momentum uncer- 
tainties that are important in the conceptual understanding of the interference phenomena. 
A wave packet description should also be implemented in the measurement or stopping (dis- 
entanglement) of the charged lepton, which we treated as an event sharp in space time at a 
time scale t c and distance L c , a wave packet treatment would smear these scales over a lo- 
calization length scale of the wave packet, which is determined by the measurement process 
(or perhaps the mean free path of the charged lepton in the stopping material). 
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FIG. 15: Same as fig. (03J with T K L C = 100. 



The typical analysis of neutrino oscillations in terms of (Gaussian) wave packets |8CH87 
clarifies that neutrino wave-packets evolve semiclassically, the center moves as the front of a 
plane wave with the group velocity and is modulated by a Gaussian envelope which spreads 
via dispersion. Wave packets associated with the different mass eigenstates separate as they 
evolve with slightly different group velocities and, when their separation becomes of the 
order of or larger than the width of the wave packet, the overlap vanishes and oscillations 
are suppressed oc e~^ L l L ^) ) w here L coh ~ a El /5m 2 and a is the spatial localization scale 
of the wave packet. This suppression becomes important when L co h < L, which for 5m 2 ~ 
eV 2 , E > 30 MeV, L ~ 100 — 600 m implies a < 1 — 5 x 10 -13 m which, while much bigger 
than nuclear dimensions, is much smaller than atomic scales. If a firm assessment confirms 
that neutrino wavepackets are produced with such localization length or smaller, then this 
decoherence effect must be introduced in the oscillation probability. 

As discussed in [80] the wave packet description also features another source of decoher- 
ence in the localization term, which suppresses coherence when a > L osc ~ E/5m 2 which is 
unlikely to be relevant in short baseline accelerator experiments. A complementary inter- 
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pretation of decoherence for Tm < 5m 2 /2E in terms of wave packets is discussed in ref. 80 
if a neutrino wavepacket produced by the decay of a parent particle of width Tm is assigned 
a localization length, a ~ 1/Tm, then the condition for decoherence from the localization 
term, a ~ L osc , becomes equivalent to Tm — Sm 2 /2E which is recognized as 1Z ~ 1 in our 
discussion. Although we do not see an obvious relation between the results obtained above 
with the non-perturbative field theoretical Wigner-Weisskopf method and the interpretation 
of a wavepacket with localization length 1/Tm, our results are certainly in agreement with 
this interpretation; however, we emphasize that the analysis above also reveals another scale 
that is important for decoherence, namely L c , which is the length scale at which the charged 
lepton that is emitted along with the neutrino is observed or absorbed. As pointed out 
above there are two important dimensionless quantities that determine decoherence in the 
plane wave limit: 71 and V mL c - 



VII. CONCLUSIONS AND FURTHER QUESTIONS. 

Motivated by the cosmological importance of new generations of heavier sterile neutrinos 
and recent proposal for high intensity sources, this article focuses on two different aspects 
related to the search of sterile neutrinos: 1) a proposal to search for heavy (~ MeV-range) 
sterile neutrinos by studying the production of negative helicity charged leptons in ir~,K~ 
decay at rest (or positive helicity in the decay of tt + ,K + ) as a complement to the search 
for monochromatic lines in the muon (or electron) spectrum, and 2) an assessment of the 
impact of decoherence effects from the lifetime of the parent meson and the stopping distance 
scale of the charged lepton on the experimental fits for sterile neutrinos masses, mixing and 
CP-violating angles in short baseline experiments. 

For searches of heavy sterile neutrinos from tc~, K~ decay at rest, we obtain the branching 
ratio for charged leptons to be produced with negative helicity (or positive helicity for the 
decay of n + ,K + ). This branching ratio determines the abundance of the negative helicity 
states in the production process and we suggest that a Stern-Gerlach type filter with a 
magnetic field with a gradient along the direction of the collimated charged lepton beam 
emitted back to back with the (anti) neutrinos allows to spatially separate the different 
helicity components. The ratio of abundances between the negative and positive helicity 
states is determined by the branching ratio (111.15ft . shown in fig. (pQ) (divided by |f/^| 2 ), 
which, in combination with the search for monochromatic lines allows, to extract both the 
mass and the element of the mixing matrix \Ui s \ 2 by fitting both the energy and abundance 
with the branching ratios. 

Upper bounds on the sterile-active mixing matrix elements from previous experimental 
searches allow us to estimate the upper bounds for the branching ratios for the different 
processes, these are given by 

Br--^ < 1(T 8 - 1(T 7 ; 3 MeV < m s < 33 MeV (VII.l) 
Br-- ev - s < HT 8 - 1(T 6 ; 3 MeV < m s < 135 MeV (VII.2) 

with the electron channel providing the largest window of opportunity because of the larger 
phase space. For K-(DAR), we find 

R r — < 1 0-9 _ 1 fl-6 fnr J 4 MeV ~ m * ~ 360 MeV ^ ~ channel ) f V TT X\ 

~ iU iU ror \ 4 MeV < m s < 414 MeV (e - channel) ' 1 J 
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These upper bounds estimates suggest that these searches could be implemented in the next 
generation of high intensity experiments. 

Short baseline experiments target new generation of sterile neutrinos in the mass range 
~ eV as suggested by the LSND, MiniBooNE results and reactor anomalies. In current 
accelerator experiments, (anti) neutrinos are produced from the decay of pions or kaons either 
in flight, as in MiniBooNE/SciBooNE, or at rest, as recent proposals suggest. We recognized 
two sources of decoherence that impact the interpretation of the data and experimental fits 
to extract masses, mixing and CP-violating angles: a) the width of the parent meson re- 
introduces an energy (or time) uncertainty and b) the stopping distance L c of the charged 
lepton that is produced in a quantum entangled state with the (anti) neutrino, decoherence 
effects are encoded in two different dimensionless quantities, 

Stti 2 - 

Tl l3 {E) = —^ ■ T M L C . (VII.4) 

The usual quantum mechanical formula for the oscillation probabilities are modified as 
follows: 



e l ^E~ — y e l - 



l — e -r M (p)ic 



1 - ifcij 



(VII.5) 



We study the impact of the decoherence effects both for Dirac and Majorana neutrinos, 
addressing in particular CP- violating effects as well as v — > V oscillations and |AL| = 2 
transitions in the case of Majorana neutrinos. In all cases, we find that, for IZij, TmL c > 1, 
the oscillation probabilities are suppressed and the oscillatory functions feature energy- 
dependent phase-shifts that results in an overall off-set that impacts the determination of 
the mass. If these decoherence effects are neglected in the experimental analysis and the 
data are fit with the usual quantum mechanical oscillation probabilities the masses, mixing 
and CP-violating angles are underestimated. 

In particular, on MiniBooNE/SciBooNE, for example, neutrinos are produced from pion 
decay for which we find 1Z ~ l/3(5m 2 /eV 2 ) and TmL c ~ 1, with one sterile neutrino with 
m s ~ 3eV, fitting with two-generation mixing underestimates sin 2 (2$) and 5m 2 by nearly 
15%. Similar underestimates follow for CP-violating angles and | AL\ = 2 processes in 3 + 2 
schemes. 

We also conclude that reactor and (current) accelerator experiments are fundamentally 
different in that the lifetime of the decaying parent particles in reactor experiments is hun- 
dreds to thousands of years, compared to pion or kaon lifetimes, and charged leptons (muons) 
are stopped within the core so that for reactors TL C , A^L C <C 1 and decoherence effects 
are all but negligible, unlike the situation for example for MiniBooNE/SciBooNE. We also 
suggest that next generation of high intensity experiments in which (anti) neutrinos are pro- 
duced from tt,K (DAR), decoherence effects may be suppressed considerably by designing 
the experiment so that the charged leptons produced with the neutrinos (mainly muons) are 
stopped on distances much smaller than the decay length of the mesons, in which case the 
usual quantum mechanical oscillation probabilities furnish an accurate description of mixing 
and oscillations. 
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Appendix A: Quantization: Mesons, Dirac and Majorana neutrinos 



We quantize the (pseudo) scalar and fermion fields in a quantization volume V. The 
charged (complex) (pseudo) scalar field is as usual 



M(x,t) = J2 



where E^f = >/p 2 + m 2 M with m M the mass of the corresponding meson. It follows that 



(A.l) 



" K *f y/2VE M (p) L P J 

It proves convenient to introduce the combinations 

M + (p,t) ± Af-(p,t) = (Ape- iEM ^^ ± (sL f e iJ5M(p)t ) . 
For Fermi fields we work in the chiral representation, 



(A.2) 
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-ex* 
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-1 



(A.3) 



(A.4) 



and for a generic Fermion /, either for charged lepton or Dirac neutrinos of mass m,f, we 
write 



(A.5) 



For Dzrac fermions of mass mj 



ip(k, h, t) 



h th U h (k) e- iE ^ + d)_ %h V h {k) e lE ^ 



(A.6) 



where Ef(k) = ^Jk 2 + m 2 and the spinors U h , Vh are eigenstates of helicity with eigenvalue 
h — ±1, these are given by 



where 



JV) = ^E f (k) + k ; e(fc) 
and v±(/c) are helicity eigenstates Weyl spinors: 



rrii 



E f (k) + k 



(A.7) 
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where 

k = k ( sin 9 cos 0, sin 9 sin 0, cos 9) . 

A useful representation is 

(1 + a-k) ( 1 



«+(*) 



V / 2(l + 



cos( 







u_(Jfe) 



(1-a-k) 



COS I 



1 



(All) 



(A.12) 



The Weyl spinors (1A.10I) satisfy 

vl(k)-v h ,(k) = V (A.13) 

Majorana fields are charge self-conjugate and generally obey 

^P c = i 7 2 iP* = e^ip (A.14) 

with £ an arbitrary (real) phase, which we choose £ = 0. In the chiral representation flA.4j) 
writing 

4>R 



if. 



it follows that 



ia 2 ip* L 
-ia 2 i>* R 



(A.15) 



(A.16) 



Therefore, a Majorana field is obtained by combining the positive frequency component with 
its charge conjugate as the negative frequency, namely 



Xfrt) EE ^2VE f {k) 



h=± 



b %h U h {k) e -Wk)t-U.S) + b\ h U£(k) e W*)*-*-*) 



where 



u:(k) = N f 



W_{k) = N f 



v + (k)^ 
e(k) v + (k) 



e{k)v-{k) 

and we have used the property 

(ia 2 )v* + {k) = -v^{k) ; (ia 2 )v*_{k) = v+(k) . 
In particular the negative chirality component of the Majorana neutrino is 



(A.17) 

(A.18) 
(A.19) 



Xl{x, t) 



k L 



E f (k) + k 
2E f {k) 



- b n + e(k) v+(k) + b L _ u_(jfe) ] e -W*) 



+ St v_{k) + St e(Jfe) v + {k) e <E f (k)t-k-Z) 



From the representation (1A.12I) . it follows that 

vl(-k) ■ v h (k) = ; h = ±. 

It is straightforward to confirm that the Hamiltonian for the Majorana fields 
1 



d 3 x X j (x,t) -ia-W + f3m f X (x,t) = J2 E f( k ) b l h k,h 



(A.20) 

(A.21) 
(A.22) 



k,h 



where the zero point energy has been subtracted. 



34 



Appendix B: Wigner-Weisskopf method for M — > Iv: 

The purpose of this appendix is to provide technical details of the Wigner-Weisskopf 
approximation as applied to the M — >• Iv process. For a more extended discussion see 
refs.@,[77l- 

The total Hamiltonian is given by H = H + Hi, where H is the free Hamiltonian and 
Hi is the interaction part. The time evolution of a state in the interaction picture is given 
by 

tj t \V(x,t)) I = H I \y(x,t)) I (B.l) 
where Hj(x,t) = e zHot Hi(x,t)e~ tHot . The formal solution of (IB . 1 1) is given by 

\*(£,t))i = U(t,t )\#(£,t )) I (B.2) 

where U(t,t ) = T(e~ 1 ^ . Expanding the state \^f(x,t))j in the basis of eigenstates 

of H we have 

|*(f,t)>/ = X) C '»(*)l n > ( B - 3 ) 

n 

where H \n) = E n \n). It is straightforward to show that J2 n \Cn(t)\ 2 = const which is a 
consequence of unitary time evolution. 

Now consider the initial state at time t = to be one meson state of definite momentum, 
namely 

|*(£,t = 0)) 7 = \M) P = J2 C ^ = ( R4 ) 



which gives the initial condition C n (t = 0) = 5. 



From eq.( 1B.ip . upon expanding in basis states, it follows that 

jC n {t) = -iJ2(n\Hi(t)\i)C m (t) (B.5) 

m 

The interaction Hamiltonian ( 111.11) connects the initial meson state, \Mp) to lep- 
tonic/neutrino states, {|/) g) These states in turn are coupled back to \Mp) via 

Hi, but also to other multiparticle states which describe processes that are higher order in 
perturbation theory. However, we will only be considering states connected to \Mp) via first 
order in perturbation theory. The case that will be of interest to us will be M — > Iv and is 
shown in Figure (fT6|) . 
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Considering the set of equations for these states, we obtain 

d 
dl 



^-C M {t) = -iY.WHm^C.it) (B.6) 



±C K {t) = -i(n\Hj{t)\M)C M {t) (B.7) 
dt 

where \n) is the intermediate state, \l a (k, h))\v a (q, h')). Using the initial conditions for t = 0, 
one obtains 

C K (t) = -i f dt\K\Hj{t')\M)C M {t') , (B.8) 
Jo 

which when inserted into ( IB. 61) leads to 

±C M (t) = - [ dt , y;<M|ff / (t)|K><«| J H J (t')|M)C w (t / ) = - I dt'H M {t-t')C M {t') (B.9) 
dt Jo K Jo 

Where the meson self energy has been introduced 

E M {t-t') = ^(Mltf^l/^liJ^OlM) = ^\(M\H I (0)\K)\ 2 S E "- E *X t -V (B.10) 

K K 

This self-energy is recognized as the one- loop retarded self energy with the \l) \u) intermediate 
state. 

Solving eq.( lB.9l) produces a solution for the time evolution of the meson amplitude. We 
can use the solution for Cm(^) to obtain an expression for the amplitudes C K (t) which allows 
for computation of the probability of occupying a particular state at any given time. We may 
solve eq.( flB.9p ) either via Laplace transform, or in the case of weak coupling, a derivative 
expansion which yields the same result at long times (t ^> 1/uim)- Here, we follow the latter 
method which is the original Wigner-Weisskopf approximation. 

We begin by defining the quantity 

W (t,t')= f dt"Z M (t-t") (B.ll) 

so that 

-^W (t,t') = E M (t-t') , W o (f,0)=0 (B.12) 
Integrating eq. (lB.9p by parts yields 

jC M (t) = - J*dt'Z M (t-t')C M (t') = -W (t,t)C M (t) + £ dt'W (t,t')-^C M (t') (B.13) 
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The first is term second order in Hj whereas the second term is of fourth order in Hj 
and will be neglected. This approximation is equivalent to the Dyson resummation of the 
one- loop self energy diagrams. Thus to leading order, eq. (1B.9j) becomes 



d 
dt 



c M (t) + w (t,t)c M (t) = o 



where 



W (t,t)= [ dt'X M (t-t')= [ 
Jo Jo 



dt'J2\( M \Hi{0)\^)\ 2 e iiE 



M -E K )(t-t') 



(B.14) 



(B.15) 



Inserting a convergence factor and taking the limit t — > oo consistently with the Wigner- 
Weisskopf approximation, we find 5 



e^0+ E M - E K + 16 2 



where 



AE M = V 



|(M|ff,(0)l")l ; 
Em — E K 



(B.16) 



(B.17) 



is the second order shift in the energy which will be absorbed into a renormalized meson 
energy and 

T M ee 27tJ2\( m \Hi(0)\k)\ 2 8(E m - E K ) (B.18) 

K 

is the decay width as per Fermi's Golden rule. Therefore in this approximation, we arrive 
at 



C M {t) = e 



-iAE M t e -^ft _ 



(B.19) 



Inserting this result into eq. (1B.8|) leads to 

C K (t) = -i(«|ffj(0)|Af) f dt'e- i{ - EM+AEM - EK - iE ^ )t ' 

Jo 



■(k\Hj(0)\M) 



1-e 



i(E M +AE M -E K -i-!f)t 



E A + AE M ~E K -l^ 



(B.20) 



Defining the renormalized energy of the single particle meson state as E r M = Em + AEm 
and passing to the Schroedinger picture \M(t))s = e~ iHot \M(t))j, we find that 



\Mz{t)) s = e 



-iHnt 



C M (t)\M) + X)C«(*)I«> 



- t e-- t |Mr(0))-^M^(0)|Mr) 



1 _ e -i(E^ M -E K -iL^)t 
^M ' 2 



(B.21) 



-iE K t 



5 The long time limit in the Wigner-Weisskopf approximation is equivalent to the Breit-Wigner approxi- 
mation of a resonant propagator 77 . 
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The interaction Hamiltonian for M — > l a u a is given by eqn. (1II.5j) and the quantization 
from Appendix A leads to the matrix element 



Ka ' A M p/ v^Y V8E M (p)E a (k)E 3 (q) 



k = p + q 



(B.22) 



which yields our final result for the entangled quantum state resulting from meson decay 



P ^ 3 ^VE M {p)E a {k)E 3 {q) 



x 



1 _ e -H^k 



(El 1 (p)-E a (k)-E J (q)-i'-f-)t 



e m(p) ~ E a (k) - Ej(q) - i^f- 



q,aj,h,h' 

-i{ Ea ( k ) +E] mtii-^k^iujih',-^)) 



(B.23) 



Appendix C: On the normalization pil.6|) : 



The normalization of the disentangled neutrino state 6 j) has another important inter- 
pretation, it is recognized as the number density of charged leptons produced from the decay 
of the meson. To see this, consider the expansion of the Dirac field for the charged lepton 
as in eqn. (IA.6[) where &t ft creates a charged lepton l~ with momentum k and helicity hi. 

The number operator for particles is = bt^b^ h , and its expectation value in the full 
meson state is given by 



n i h =^mn k * M \Mz{t)))= y: 



U aj\ 2 \M^j{k,q, hi,h')\ 2 



f 8VE M (pEj (q)E a (k) 



Tajik, q;t c ] (q;h 



which is recognized as the normalization (1III.6I) . namely 

K(g; hi) = 4 A . 



(C.l) 



(C.2) 



From the definition of the partial width ^m-^i^vj(Pj ^i, h!) of meson decay into a lepton 
a of helicity h and neutrino eigenstate z/j of helicity h' 



1 



d 3 q \MZAk,q,hi,h')\ 
9 ' ° A ,l 27r5(E M ( J) )-E a (\p-q\)-E j ( q )) 



2E M (p)J (2tt) 3 2E j (q)2E a (k) 
and the total decay width 

r A/(p) = \ U °3 ] ? V M--4GV j fa K h') 
j,h' 



(C.3) 
(C.4) 



it follows that the total number of charged leptons produced at time t c is given by 

d 3 q 

h J (2^p ' " 1 ^ I i 



/ f>-)3 k,hi J 



N v (q,hi 



I _ e ~ r M(p)tc 



(C.5) 
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